Physical Bounds to the Entropy-Depolarization Relation in Random Light Scattering 
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We present a theoretical study of multi-mode scattering of light by optically random media, 
using the Mueller-Stokes formalism which permits to encode all the polarization properties of the 
scattering medium in a real 4x4 matrix. From this matrix two relevant parameters can be extracted: 
the depolarizing power Dm and the polarization entropy Em of the scattering medium. By studying 
the relation between Em and Dm, we find that all scattering media must satisfy some universal 
constraints. These constraints apply to both classical and quantum scattering processes. The 
results obtained here may be especially relevant for quantum communication applications, where 
depolarization is synonymous with decoherence. 
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Introduction The polarization aspects of random 
light scattering have drawn quite some interest in re- 
cent years, since they present a diagnostic method of the 
medium involved and also help visualization of objects 
that are hidden inside the medium When polarized 
light is incident on an optically random medium it suf- 
fers multiple scattering and, as a result, it may emerge 
partly or completely depolarized. The amount of depo- 
larization can be quantified by calculating either the en- 
tropy (Ef) or the degree of polarization (Pf) of the scat- 
tered field Q. It is simple to show that the field quanti- 
ties Ef and Pp are related by a single-valued function: 
Ef = Ep(Pp). For example polarized light (Pf = 1) 
has Ef = while partially polarized light (0 < Pp < 1) 
has 1 > Ef > 0. When the incident beam is polarized 
and the output beam is partially polarized, the medium 
is said to be depolarizing. An average measure of the de- 
polarizing power of the medium is given by the so called 
depolarization index (Dm) S3J- Non-depolarizing media 
are characterized by Dm = 1, while depolarizing media 
have < Dm < 1- A depolarizing scattering process is 
always accompanied by an increase of the entropy of the 
light, the increase being due to the interaction of the field 
with the medium. An average measure of the entropy 
that a given random medium can add to the entropy of 
the incident light beam, is given by the polarization en- 
tropy Em H- Non-depolarizing media are characterized 
by Em — 0, while for depolarizing media < Em < 1. 
As the field quantities Ep and Pp are related to each 
other, so are the medium quantities Em and Dm with 
the key difference that, as we shall show later, Em is a 
multi- valued function of Dm- 

The purpose of this Letter is to point out a universal 
relation between the polarization entropy Em and the 
depolarization index Dm valid for any random scatter- 
ing medium. This relation covers the complete regime 
from zero to total depolarization. It has been intro- 
duced before, by Le Roy-Brehonnet and Le Jeune Q, 
in an empirical sense, to classify depolarization mea- 
surements on rough surfaces (sand, rusty steel, polished 



steel, . . . ) . We derive here its theoretical foundation and 
present analytical expressions for the multi-valued func- 
tion Em = Em(Dm)- Although the (Em, Dm) relation 
is essentially classical, we use a single-photon theoreti- 
cal approach, exploiting the well known analogy between 
single-photon and classical optics |f| . We prefer this to 
a classical formulation since it offers a natural starting 
point for extension to entangled twin-photon light scat- 
tering by a random medium, which is a true quantum 
phenomenon that could deteriorate quantum communi- 
cation. 

Polarization description of the field Let us consider a 
collimated light beam propagating in the direction z. In a 
given spatial point r, the monochromatic time-dependent 
electric field associated with the beam is a complex- 
valued vector E(i) = X(t)x + Y(t)y. This vector de- 
fines the instantaneous polarization of the light which 
is, in any short enough time interval, fully polarized. Al- 
ternatively, the same light beam may be described by 
a time dependent real- valued unit Stokes vector s(t) = 
{2Re(X*Y),2lm(X*Y), \X\ 2 -\Y\ 2 }/(\X\ 2 + \Y\ 2 ), which 
moves on the Poincare sphere (PS) |fj. Of course, no de- 
tector can measure the instantaneous polarization, the 
best one can get is an average polarization over some 
time interval T. If during the measurement time T the 
Stokes vector s(t) maintains the same direction, then the 
beam is polarized. Vice versa, if s(t) moves over the PS 
covering some finite area, then the beam is partially po- 
larized. In the last case, for stationary beams, the mo- 
tion of s(t) produces a probability distribution over the 
PS which determines the degree of polarization of the 
light • Time dependence of the polarization is not the 
only cause for depolarization, also spatial dependence, 
for example, may lead to loss of polarization. 

We stress that this picture is not limited to the clas- 
sical domain; in Ref. || we found, e.g., that a multi- 
mode single-photon scattering process generates a k- 
dependent Stokes vector distribution. More generally, if 
ip = {t, k, A, . . .} denotes the set of all variables (e.g., time 
t, momentum k, polarization A, . . . ) on which s — s(ip) 
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depends, then the state of a polarized light beam (either 
classical or quantum), may be described by a 2 x 2 matrix 
p(ip) = (<To + s(ip) ■ c)/2, where a is the 2x2 identity 
matrix and er = {a±, ct 2 , 03} arc the Pauli matrices. The 
matrix php) is known as the coherency matrix in classical 
optics |(j and as the density matrix in quantum mechan- 
ics 0. Since by construction Trp(^) = 1, each matrix 
p(ip) can describe either a purely polarized beam in clas- 
sical optics, or a pure photon state in quantum optics 
[lf1 |. However, the state of a partially polarized beam 
must be described by the matrix p — J [dtp) w(ip)p(ip), 
where J (dtp) is the integration measure [Tl( in the space 
of the variables ip and J(dip)w(ip) = 1- The statistical 
weight w(ip) > 0, defines a probability distribution over 
the PS. It is clear that p can represent a mixed photon 
state in the context of quantum optics as well. If A de- 
notes any polarization-dependent observable, its average 
value must be calculated as: 



matrix ||Mu^|| appears which is defined as 



(A)=Tr(pA) = / (#MVWMV0). 



(1) 



If A represents the entropy of the field, i.e. A = — log(p), 
then (A) = — Tr(plogp), which is the von Neumann 
entropy S of the photon state However, by us- 

ing Eq. Q it is easy to see that this coincides with 
the Gibbs entropy 13] of the distribution w(ip), since 
S = — J (dip) w(ip) log(w(ip)), in agreement with the re- 
sults of Ref. 0. 

Single-photon scattering and multi-mode Mueller for- 
malism The theoretical framework for studying one- 
photon scattering has been established elsewhere [8| , here 
we use the results found in |£| to extend the Mueller- 
Stokes formalism to quantum scattering processes. In 
classical optics a polarization scattering process can be 
characterized by a real-valued 4x4 matrix, the so 
called Mueller matrix M [2|, which describes the po- 
larization properties of the scattering medium. Wc 
show now that such a matrix description can be ex- 
tended to the quantum (single-photon) scattering case. 
Let us consider a photon prepared in the pure state 
p(ip), approximatively described by a monochromatic 
plane wave |ko,Ao). In this case ip = {ko,Ao}. Now, 
let us suppose that the photon is transmitted through 
a linear optical system described by an unitary scat- 
tering operator T such that p(ip') = p(ip)T repre- 
sents the pure state of the photon after the scatter- 
ing, where ip' is the set of all scattered modes: ip' = 
{ki, Ai, k2, A2, • • ■}■ A multi-mode detection scheme im- 
plies a reduction from the set ip' to the subset of the 
detected modes ip" — {ki, Ai, . . . , kjy, Ajv} C ip' which 
causes a transition from the pure state p(ip') to the mixed 
state p — J (dtp") w(ip")p(ip"). If we denote the Stokes 
parameters of the beam before and after the scattering 
with = Tr(p(ip)a fl ) and s' respectively (p = 0, 1, 2, 3), 



then the classical result s' 
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M^ v s v is retrieved, 



M^ v oc / dkm A1 „(k). 
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The local (with respect to the momentum) matrix ele- 
ments m^v (k) are defined by means of the matrix relation 

^ T (k)T(k,koKTt(k ,k)^(k) = m pv (k)a u , (3) 

(p, v — 0,1,2,3), and summation over repeated indices 
is understood. Explicit expressions for the 2x2 matrices 
W(k) and T(k, k') can be found in Ref. |||. The propor- 
tionality factor in Eq. <0 can be fixed by imposing the 
condition Moo = 1- When ip" reduces to a single mode 
{k, A}, then Wij(k) = Sij and the classical formalism is 
fully recovered. 

Depolarization index Dm and polarization entropy Em 
Now that we have a recipe to calculate the Mueller 
matrix describing a multi-mode scattering process, we 
use this knowledge to study the depolarization prop- 
erties of the scattering medium. Within the Mueller- 
Stokes formalism, the degree of polarization Pp of the 
field and the depolarization index Dm of the medium, 
are defined as P F = (sf + s\ + s^) 1 / 2 / s and D M = 

(Tr(M T M)/3 - 1/3) 1/2 , respectively, where s M [p = 
0,1,2,3) are the Stokes parameters of the field and 
Moo = 1 has been assumed. A deeper characterization 
of the scattering medium can be achieved by using the 
Hermitian matrix H [pH fl5| defined as 
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M M „ (07, ® cr*) , 



(4) 



where Tr(H) = 1. It can be shown |4| that a physically 
realizable optical system is characterized by a positive- 
semidefinite matrix H . Let < A„ < 1 (v = 0, . . . , 3) 
be the eigenvalues of H. Then it is possible to express 
both the depolarization index Dm and the polarization 
entropy Em as a function of the A„'s. Explicitly we have 
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1/2 



and 



E M = A„log 4 (A„). 



(5) 



(G) 



v=0 



with the difference that a generalized (measured) Mueller 



Now we are ready to show the universal character of 
the (Em, Dm) plot originally introduced in Ref. 0]. 
More precisely, we show that it allows to characterize 
all possible scattering media by means of their polari- 
metric properties. The main idea is the following: both 
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Em and Dm depend on the four real eigenvalues of H 
which actually reduces to three independent variables be- 
cause of the trace constraint Tr(i?) = 1. If we use Eq. 
(J3J) to eliminate one of these variables in favour of Dm 
we can write Em = Em{Dm,ol, (3) where a, (3 represent 
the last two independent variables. Then, for each value 
of < Dm < 1, different values of Em can be obtained 
by varying a and (3 between and 1. In such a way we 
obtain a whole domain in the Dm-Em plane instead of 
just a curve. In order to do that, we have implemented 
a Monte Carlo code to generate a uniform distribution 
of points over the 4-dimensional unit sphere: the square 
of the four coordinates of each point is an admissible set 
of eigenvalues of H . In this way we have generated the 
graph shown in Fig. 1. The boundary of this domain 



1 

0.8 
0.6 



^1 



0.4 
0.2 








; p 2 












0.2 


0.4 0.6 0.8 1 




D M 



FIG. 1: Numerically determined domain in the Dm — Em 
plane corresponding to all physically realizable polarization 
scattering processes. The solid curves are the analytically 
obtained bounds. The four cusp points pi — (0, 1), p2 = 
(1/3, log 4 3), ps — [l/y/i, 1/2), pi — (1,0) separate different 
polarization scattering processes, as described in the text. 

is formed by the curves — 1, ... ,4), joining the 

points (pi — * pj). The analytical expressions for these 
curves are 

E(n, /) = - [(1 - „/) log 4 (1 - nf) + nf log 4 (/)] , (7) 
where 

1 



n + 1 



3n + l, o 
1 ± i/l - -— — (1 - D 2 



4 n 



Ml 



(8) 



The links between the functions E(n, /) and the curves 
Cij are given in Table I where we have defined E13 = 
— (1 — /i) log 4 (-!-^ ± ) — julog 4 (^). The curve C14 is special 
in the sense that it sets an upper bound for the entropy 
of any scattering medium. We find numerically that the 



TABLE I: List of the analytical curves (continuous lines) in 
Fig. 1. The second column refers to the equations generating 
the corresponding curves, while the third column gives the 
eigenvalues of H. The first four curves form the boundary of 
the physical domain; the last two represent inside curves. 



Curve 


Generating equation 


Eigenvalues of H 


Cia 


E(3,f+) 


{X,u,^,u} 


C23 


E(2,f+) 


{\,H,H,0} 


C34 


E(hf±) 


{X,u,0,0} 


C14 


E(3,f-) 


{A, /j,, p,, fj} 


Cl3 


E13 


{X,X,u,u} 


C24 


E(2,f-) 


{\fj,,fi, 0} 



value of the entropy on this curve is very well approxi- 
mated by 



E 



M 



1 



U M) ■ 



(9) 



where 7 = 0.862, which is, interestingly, almost equal to 
e/n. Then, for all depolarizing scattering media the con- 
dition E M < E M must be satisfied. It is interesting to 
note that a purely depolarizing scattering medium (with 
diagonal Mueller matrix) leads to Em = Ec- By using 
thermodynamics language, one may interpret Fig. 1 as 
a polarization "state diagram" where different phases of 
a generic scattering medium, characterized by different 
symmetries of the corresponding Mueller matrices, are 
separated by the curves Cij. It is worth to note again 
that there is nothing inherently quantum in the above 
derivation of the physical bounds Eq. JJJ , therefore these 
results have validity both in the classical and in the quan- 
tum regime. 

Random matrix approach We have checked the va- 
lidity of the theory outlined above, for scattering media 
in the regime of applicability of the random-matrix the- 
ory (RMT) 0| . Random media, either disordered me- 
dia [l| or chaotic optical cavities jljj , can be represented 
by ensembles ^^|. The transmission of polarized light 
through a random medium may decrease the degree of 
polarization in a way that depends on the number N 
of the detected modes via Eq. J5J). Under certain con- 
ditions, RMT can account for a statistical description 
of the light scattering by random media 0, |2(|. Let 
be the complex probability amplitude that a pho- 
ton is scattered in the state |k, A). Then, according to 
RMT, the real and the imaginary parts of the scatter- 
ing amplitudes i/>A(k) are independent Gaussian random 
variables with zero mean and variance that can be fixed 
to 1. The assumption of independent variables is jus- 
tified since usually the set ip" of the detected modes is 
much smaller than the set ip' of the all scattered modes 
|2l| . Let us suppose now that the impinging photon is in 
the pure state |k , Ao). In this case ifrxQc) — T\\ (k, k ) 
and the statistical distribution of the 's can be nu- 
merically calculated accordingly to Eqs. (l^l.'St . In this 
way we have calculated the ensemble-averaged polariza- 
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tion entropy (Em) and depolarization index (Dm) of the 
medium, as functions of N for the case in which the angu- 
lar aperture of the detector is so small that (k) ~ S{j . 
The results are shown in Fig. 2 for the cases of a generic 
scattering medium (73a (k, ko) unconstrained) and of 
a polarization-conserving medium (T\\ (ls., ko) oc S\\ ). 
The last case is realized when the geometry of the scat- 
tering process is confined in a plane. As one can see, 
for both cases RMT results cover only a small part of 
the (Em, Dm) diagram; however, the numerical data are 
consistent with the analytical bounds given by Eq. Q. 
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FIG. 2: RMT results for the ensemble-averaged polar- 
ization entropy (Em) as function of the ensemble-averaged 
depolarization-index (Dm) for generic (dark squares) and 
polarization-conserving (open squares) scattering processes. 
In both cases each point correspond to a given number N of 
detected modes. When N increases from 1 to 30, points move 
from the bottom to the top of the figure. The solid lines are 
the analytical bounds of Fig. 1. 

Conclusions In summary, we have studied the scat- 
tering of light by optically random media, from a polar- 
ization point of view. To this end we have first extended 
the Mueller-Stokes formalism to make it suitable for the 
description of single-photon scattering processes. Then, 
after the calculation of the Mueller matrix M characteriz- 
ing the polarization properties of the scattering medium, 
we have extracted from M the depolarization index Dm 
and the polarization entropy Em- By analyzing the func- 
tional relation between Em and Dm , we have found that 
the polarization properties of any scattering medium are 
constrained by some physical bounds. These bounds have 
an universal character and they hold in both the classical 
and the quantum regime. Our results provide a deeper 
insight in the nature of random light scattering by giving 
an useful tool, both to theoreticians and experimenters, 
to classify scattering media according to their polariza- 
tion properties. The use of this tool may be particularly 



relevant in quantum communication where it is desirable 
to manipulate and control the polarization of the light 
|22|. Presently, experiments are in progress in our group 
to verify this theoretical framework. 
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